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We extend the well-known Cont-Bouchaud model to include a hier¬ 
archical topology of agent’s interactions. The influence of hierarchy on 
system dynamics is investigated by two models. The first one is based on 
a multi-level, nested Erdos-Renyi random graph and individual decisions 
by agents according to Potts dynamics. This approach does not lead to 
a broad return distribution outside a parameter regime close to the orig¬ 
inal Cont-Bouchaud model. In the second model we introduce a limited 
hierarchical Erdos-Renyi graph, where merging of clusters at a level h + 1 
involves only clusters that have merged at the previous level h and we use 
the original Cont-Bouchaud agent dynamics on resulting clusters. The sec¬ 
ond model leads to a heavy-tail distribution of cluster sizes and relative 
price changes in a wide range of connection densities, not only close to the 
percolation threshold. 

PACS numbers: 89.65.-s, 89.65.Gh, 89.75.-k, 05.40.-a 


1. Introduction 

The paradigm of investor’s rationality and effectiveness was ousted in 
the second half of the 20th century by the behavioural economics. A mile¬ 
stone was the work of Kahneman and Tversky written in 1979 [T]. Since 
then, many studies were done and many models describing the stock market 
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were developed [SlEliaElElElElE]. One of the first microscopic model 
was proposed by R. Cont and J.P. Bouchaud [TO]. It was suggested that 
the heavy tails observed in distributions of stock returns, corresponding to 
large fluctuations in prices, arise as a result of collective phenomena such 
as herding behaviour. In the Cont-Bouchaud (CB) model, the groups (clus¬ 
ters) of investors make a collective decision during each time period and 
may choose either to buy the stock with the probability P{(p = -|-1) = a, to 
sell it with the same probability P{(l) = —1) = a, or to stay inactive with 
the probability P{4> = 0) = 1 — 2a. The clusters are formed according to the 
theory of Erdos-Renyi (E-R) random graphs: for any pair of nodes (agents), 
let p = c/N be the probability that these nodes are linked together, where 
N is total number of nodes and c is a positive parameter, which represents 
the willingness of agents to form groups. The demand/supply created by a 
cluster of investors is proportional to its size. The price of an asset changes 
from one time step to another and the relative price change is proportional 
to the total excess demand. Since for parameter c close to and smaller 
than 1, the cluster size distribution decreases asymptotically as a power law 
with exponential cutoff, the model can be solved analytically in the limit 
2aN ~ 1 according to the percolation theory. Many modifications and ex¬ 
tensions of the original CB model were proposed in past years: changing 
the topology from the random graph to a square {d = 2), cubic {d = 3) or 
hypercubic lattice (d = 4 7) [H], introducing moving agents [12] or fun¬ 

damentalists m, introducing the dependence between the activity a and 
the cluster’s size m or between a and a current price m and much more 
[ISKIZKIBlliaEniEIKSS]. 

All above studies assume only one level of the intimacy between traders 
(they are connected or they are not) and an infinitely strong interactions 
inside clusters. These assumptions are hidden in oversimplified topologies 
like the random graph or the square lattice. Our aim was to adapt the 
main idea of the CB model to hierarchical topologies. We expect that com¬ 
munity of market investors, like other social systems, exhibit hierarchical 
structures. An investor can divide others into different levels of ’’closeness” 
- small group of friends or colleagues may influence their decisions in major 
way, including agreeing on common strategy, while the direct influence of the 
community as a whole may be small. There is rarely collective behavior of a 
whole market, while this may happen much more often in close-knit groups. 
Our model is a generalization of Cont-Bouchaud approach that considered 
only binary ’’deciding together” or ’’deciding separately” interactions be¬ 
tween agents. The nature of the interactions may be also interpreted as 
access to a common source of information, which translates into follow¬ 
ing the same trading strategy in an attempt to exploit that information. 
Considering such interpretation one may argue that our models describe 
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inefficient market, where dynamics can be sometimes dominated by large 
groups following single information sources, meaning it does not follow the 
aggregate of all information or may be biased. Topological features of hi¬ 
erarchical systems have been observed in real networks [23], and these kind 
of systems have been considered in several other studies, including social 
dynamics [21] and information propagation [25]. We developed two models, 
which use concept of hierarchy for extending the original CB model in two 
different ways. 


2. Hierarchical Erdos-Renyi graph 

Consider a topology which consists of N nodes which are grouped into 
many nested clusters. The procedure of creating such a complex system is 
following: 

1. Link together each pair of nodes with a probability pi = c/N. The 
result is Wi clusters with the hierarchical degree h = 1. The nodes 
belonging to the same cluster are nearest neighbours, so their mutual 
degree of neighbourhood is one. 

2. The existing clusters of hierarchical degree h = \ are treated as 
primary-level nodes which are linked together randomly with a prob¬ 
ability p 2 = c/Wi per pair. The result is W 2 clusters with the hier¬ 
archical degree h = 2, containing sets of clusters of degree h = 1. If 
nodes i and j belong to the same cluster with the hierarchical degree 
h = 2, but to different clusters of the hierarchical degree h = 1, their 
mutual degree of neighbourhood is two. 

3. Repeat step 2 for following hierarchy levels h, until hierarchy level 
h = H 'is reached, where all the nodes belong to only a few clusters. 
The number of iterations H required is proportional to logarithm of 
system size H ~ In A^. Assume all the clusters of level H form a single 
cluster of hierarchical degree H + 1 that contains all nodes. 

The resulting graph will be called a hierarchical Erdos-Renyi graph (HERG). 


3. Cont-Bouchaud model on hierarchical Erdos-Renyi graph with 

Potts interactions 

We define a hierarchical Cont-Bouchaud model (HCB) in the following 
way. Let nodes of HERG represent the stock traders. At each time step, 
each agent has to take a decision about trading: buy {(j) = 1), sell {(j) = — 1) 
or do nothing ((/> = 0), just as in the standard CB [TU] model. 
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Fig. 1. The structure of hierarchical Erdos-Renyi graph (HERG). The nodes in red 
clusters are the nearest neighbours. Green clusters represent the neighbourhood of 
the third degree. Clusters are formed through connecting clusters of previous level 
at random, as in E-R random graph. 


The agents are picked out at a random order and take decisions in¬ 
dividually, but not independently, with probabilities corresponding to the 
probabilities of spin directions in Potts model in the canonical ensemble. 
The probabilities are 

P{ct>i) = ( 1 ) 

where Z = -p is the partition function and 

(3 = l/ksT. E{(f)i) is the spin energy corresponding to Potts Hamiltonian 
with interaction strengths Jij dependent on mutual degree of neighborhood 
h{i,j) of interacting spins: 


E{4>i) = - ^ - B6{(j)i,0), (2) 

where 6{(pi,(j)j) is the Kronecker delta, and the sum is over all agents j ^ i. 
The parameter B acts as an external held which restrains (or enhances if 
it is negative) trading and therefore plays the same role as the parameter a 
in CB model. In fact both parameters are related as a = l/(e^^ -|- 2), so 
that when an agent is not inhuenced by others, then P(-|-l) = P(—1) = a. 
The coupling factor Jh = Ja^~^ depends on degree of neighbourhood h 
and J = Ji is the coupling constant which is equal to the coupling factor 
of the nearest neighbours. The parameter a < 1 controls how quickly the 
interactions weaken with the degree of neighbourhood h. At the beginning 
of each time step all agents are in a transient state null, and do not influence 
other agent decisions. During each time period the agents are appointed to 
make decisions in a random order, interacting according to EqU but only 
with agents in non-null state that already made their decisions in this time 
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Fig. 2. Tails of return distribution densities for HCB model for N = 5000, c = 0.99, 
a = 10“®, /3 = 1, J = 1000 and B pe 4.5 (corresponding to a = 0.01 in CB model). 
Distributions in BCB model for this parameter set are the same as for the CB 
model however deviations from that parameters result in fat tails disappearing. 
Top left: comparison of BCB (black cross) with the original CB model (orange X) 
with the same parameters (a = 0.01). Top right: for decreasing field B: « 4.5 
(black cross), 3.0 (blue square), 1.5 (red circle) and 0.0 (green triangle). Bottom 
left: for increasing a: 10“® (black cross), 0.1 (blue square), 0.5 (red circle) and 
0.7 (green triangle). Bottom right: for decreasing J: 1000 (black cross), 3.0 (blue 
square), 2.0 (red circle), 1.0 (green triangle). The grey line is for the visual guidance 
only and it correspond to a power law with the exponent —2.5. 


step. Then the return is calculated in a similar way as in CB model: 

N 

Ax oc'^(j)i{t). (3) 

i=l 

Let us note that there is no memory in agent states 4>i that are reset to the 
state null after each time step. 

We have investigated the model numerically for various parameters a, J, (5 
and B, starting from the parameter set corresponding to the original CB 
model, where we assumed a = 10“®, J = 1000, B ~ 4.5 (which correspond 
to o = 0.01) with (3 = 1. 
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When parameters are chosen so that the HCB model is replicating the 
original CB model, the return distribution is the same as resulting from the 
CB dynamics. This behavior, along with results of changing parameters 
can be seen on Fig. [2l where the data have been aggregated over r = 100 
realization with t = 1000 time steps each. When the HCB model parameters 
deviate from those reflecting the CB model, the system stops working in the 
“percolation regime” where the return distribution reflects the cluster size 
distribution and instead it works in the “Potts regime” where it follows the 
fluctuations arising from Potts interactions between agents. We have not 
observed fat tails in the return distribution in the “Potts regime” for any 
analyzed parameter ranges (not only varying a single parameter from the 
CB parameter set). 

4. Limited hierarchical Erdos-Renyi graph 

In the Sect. [2] we introduced a hierarchical Erdos-R&yi graph where 
all clusters of a hierarchical degree h were considered as nodes at the next 
hierarchy level h + 1 and they could be randomly linked to form clusters 
of that hierarchical degree. Here we introduce a limited hierarchical Erdos- 
Renyi graph (LHERG). The difference as compared to HERG model is that 
the merging of clusters at the level /i + 1 involves only clusters that have 
merged at the previous level h. In effect the cluster growth is limited only 
to clusters that succeed in finding partners at previous hierarchy level. 



Fig. 3. Multilevel growth of clusters in LHERG model. Random E-R connections 
between single nodes (black) form clusters of level h = 1 (orange, violet, green). 
Random connections between those clusters (single nodes do not participate, since 
they failed to connect) form clusters of level h = 2 (red). Note that after the 
second step, single nodes and the orange cluster are excluded from the further 
merging. We forget about the internal structure of clusters when we use the re¬ 
sulting disconnected clusters (orange, red and single nodes) for CB model in Sect. 
0 


At the beginning, we consider the system of N independent nodes that 
we treat as “level zero” clusters of size 1. The procedure is following: 

1. Link together each pair of nodes with the probability pi = c/N. The 
result is Wi clusters which advance to the next step. Not linked nodes 
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are clusters of size one and they do not participate in further steps 
(and they are not included in Wi). 

2. Merge together each pair of clusters with the probability p 2 = cjWi. 
The result of two clusters merging is a new cluster whose size is the 
sum of sizes of the merged clusters. New clusters {W 2 ) advance to the 
next step. The clusters which did not merge during this step do not 
advance and stop growing (and they are not included in W 2 )- 

3. Repeat the step 2 for clusters of successive levels until all clusters stop 
growing. 
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Fig. 4. The relations between the highest level of hierarchy F[ and the size of the 
system N (left picture) and between iJ and the parameter c (right picture) for 
LHERG model. 


Denote the average number of steps of the procedure of cluster growth 
by H(c,N), which can be interpreted as the highest level of hierarchy of 
LHERG. H{c, N) is proportional to the logarithm on N and has a maximum 
as a function of c (this can be seen in Fig. 0]). Note that LHERG model 
results in many disconnected clusters (resulting from not connecting clusters 
that have been dropped out of the merging procedure), unlike in HERG 
model where all agents belong to a single top-level cluster. 


5. Cont-Bouchaud model on the limited hierarchical Erdos-Renyi 

graph 

Now we apply the dynamical rules of the standard GB model to the 
clusters obtained in the LHERG model. This model is significantly different 
from the model presented in Sect. [3l and the differences can be expressed 
in 2 points: 

1. The first model retains the hierarchical structure of the HERG clus¬ 
ters, while in the second it is forgotten after the clusters are deter¬ 
mined through LHERG growth procedure. Note that LHERG results 
in multiple disconnected clusters. 
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2. The first model uses individual agent decisions and Potts interactions, 
while the second one uses collective random decisions by whole clusters 
as in the original CB model. 




Fig. 5. Comparison between the original CB model (left picture) and the CB model 
on the LHERG clusters (right picture). Distributions of price changes are computed 
for parameter c = 0.99 (black cross), 0.9 (blue square), 0.8 (red circle), c = 0.7 
(green triangle). The grey line is for visual guidance and has exponent —2.5. In 
the original CB model distribution quickly loses fat tail when c goes away from 
percolation threshold, while in our LHERG model the power-law prevails. 

Numerical simulations show that our second model displays distributions 
of price changes (returns) with heavy tails for a wide range of parameter 
c, as seen in Fig. O This means that assuming a hierarchical organization 
according to LHERG model, it is possible for CB model to explain fat tails 
of return distributions without making strict assumptions about the con¬ 
nection mechanisms between agents, such as operating at the percolation 
threshold. We would like to note that while the distibution shape becomes 
insensitive to the parameter c, the order flow parameter a changes return 
distributions just as in the standard CB model. 

The CB model on LHERG clusters produces power-law fat tails when the 
results are averaged over multiple realizations. When we keep to one, single 
realization, the distribution does not display power-law tail, and instead 
exhibits secondary peaks, that correspond to large clusters taking decisions 
to buy or sell (Eig. [6]). This means that the model can reproduce fat tails 
only when the topology is dynamical and changes in time, although not 
necessarily every time step. 


6. Conclusions 


We have introduced a model of hierarchical Erdos-Renyi graph, where 
nested clusters of successively higher hierarchies emerge from connecting 
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Fig. 6. Distribution of price changes for CB model on single realization of LHERG 
cluster. Distribution are computed for parameter c = 0.99 (black cross), 0.95 (blue 
square), 0.8 (red circle), 0.5 (green triangle). 


lower-level components. We have studied an extension of Cont-Bouchaud 
model, where we use interactions drawn from Potts model to decide whether 
agents buy or sell, allowing for different couplings between agent decisions, 
rather than an absolute correlation assumed in the CB original model. We 
have used the hierarchical E-R graph as the topology, with couplings de¬ 
creasing with a higher hierarchical degree of neighborhood. The model does 
show a broad return distribution only when working in the limit resembling 
the original CB model, meaning that the hierarchical structure of interac¬ 
tions does not result in power-law fluctuations when Potts interactions are 
used. In the second approach the CB model has been studied on clusters 
taken from the limited hierarchical Erdds-Renyi graph, where clusters can 
merge in following hierarchy levels only if they successfully make at least 
one connection in the current one. If such a structure is kept dynamics, 
the cluster size distributions show on average power-law behavoir, therefore 
inducing a power-law return distribution in CB model using this topology. 
This happens for a wide range of average degrees c, not only when it is close 
to the E-R percolation threshold, therefore lifting one of the assumptions re¬ 
quired for the model to explain fat tails in market fluctuation distributions, 
although necessitating dynamical cluster structure. 
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